Не пугайтесь странному формату этой страницы. Дело в том, что материалы на ней такие старые, что я что-то не смог придумать, как же их более внятно оформить. Речь пойдёт об одной программе, написанной ещё в 1994-м году, и о результатах её работы.
Просто, проводя в очередной раз уборку на своей машине, я раскопал кучу программ, написанных на Паскале начиная с 91-го года, и решил попробовать их запустить. Очень интересно и полезно, кстати, читать свой код многолетней давности и изучать, как изменилось собственное мышление за такой период времени.
Turbo Pascal 7 сегодня -- freeware, так что с загрузкой и установкой его проблем никаких не возникло. Запуск, правда, был не столь гладок: CRT-шный модуль отказывался инициализироваться на слишком быстрых современных машинах, но для этой проблемы нашёлся patch, и вот, о чудо, в 2009-м году на 2003-м Windows Server-е взлетает (в окошке!) Turbo Pascal, компиллирует программы, написанные мною ещё в начале 90-х, и исполняет их!
Не правда ли, занятно? Вот Вы читаете эту страницу с моего сервера, а он в это время, быть может, Паскаль крутит? :)
Впрочем, рассказ о другом. О том, как в далёком 94-м году Денни Гурский написал программку для расчёта и отрисовки множеств Жюлиа. Очень красивые получались картинки. Ну а я занялся тем, что с её помощью исследовал поведение этого множества на ряде других аналитических функций, из чего получались ещё более интересные картинки.
Настоящим исследованием это назвать нельзя (в лучшем случае -- лишь полётом математической интуиции), да и строгости в нём нет никакой. Может быть, мы и могли бы довести его до ума 15 лет назад, но не довели, и всё, что сегодня осталось -- это идея и множество цветных картинок.
О них, собственно, и пойдёт речь.
В программе множество Жюлиа на комплексной плоскости определялось как множество таких точек Z, в которых последовательность отображений
Zn+1 = Zn2 + c
не стремится к бесконечности при некотором заданном комплексном c при n стремящемся к бесконечности. Это определение, вообще говоря, несколько отличается от классического, но у нас в те времена не было Википедии и Интернета, чтобы свериться, не вылезая из-за компутера :)
На практике, разумеется, до бесконечности никто не считал. Если последовательность ещё не расходилась при некотором "пороговом" n, точка окрашивалась в чёрный цвет. Если расходилась -- в другой, в зависимости от номера шага, на котором это произошло. Расходимость определялась как превышение |Z| некоторой константы.
(Заметьте, что всё это вычислялось на 286-х машинах, так что алгоритм вынужденно был очень простым. Ибо даже картинка 200x200 могла отрисовываться весьма небыстро. Сегодня тот же самый код работает раз в 100 шустрее).
То, что получалось, было изображением даже не самого множества Жюлия, но точек, к нему прилегающих. Тем не менее, в зависимости от выбора c результат мог выглядеть достаточно красиво и очень по-разному. Например, вот что выдавалось при c = -0.82i:
Это множество фрактально. То есть, разглядывая его со сколь угодно большим увеличением, мы будем обнаруживать элементы, во всём подобные по форме более крупным. И так вглубь до бесконечности.
Вот, скажем, центральная часть показанной выше картинки:
Если увеличить её ещё сильнее, можно разглядеть маленькие "сходящиеся мосты", подобные большому в центре (они выделены белыми квадратами):
Ещё в 93-м году меня заинтересовал вопрос: а обязательно ли использовать для построения множества Жюлиа исключительно лишь квадратичное отображение? Ведь очевидно же, что этот процесс обобщается как минимум на четвёртую степень (ибо это есть ни что иное, как то же квадратичное отображение, где добавление константы через шаг пропускается). При этом даже критерий расходимости менять не надо, он остаётся строго тем же.
И вот году в 94-м я взял и попробовал. Действительно, для Zn+1 = Zn4 + с при c = 0.82i получается вот такая картинка:
Действительно, при N = 11 получается 11 сегментов (квадрат 1.6 на 1.6, c = -0.9i):
Проверка показывает, что тоже получается:
Здесь использовано отображение Zn+1 = Zn7 + Zn2 + c при c = -0.7i.
Занятно, что 7-ми и двухлучевая симетрии как будто соревнуются за влияние и в итоге разделяют его сферы: крупным структурам достаётся 180-градусная симметрия, а мелким -- 7-ми-лучевая.
Zn+1 = exp(Zn)+ c при c = 3.1i, квадрат 10 на 10 в центре координат:
С эстетической позиции не особо впечатляет... но! Есть определённое сходство с картинкой про 11-ю степень. Это же... круг с центром в минус полубесконечности, и радиусом в половину бесконечности! И, на первый взгляд, тоже фрактальный. Ага! Попробуем покрутить константу?
Константа равна 3i.
Нисколько не удивительно: экспонента периодична вдоль мнимой оси, косинус -- вдоль действительной. Поскольку их ряды родственны, то и картинки получаются в некоторой степени похожими. Но косинусовские мне нравятся больше. Особенно если поиграть константами:
Zn+1 = Ln(Zn)+ c при c = -1.2i, квадрат 10 на 10:
Не впечатляет... :( Ну то есть повторяемость наблюдается, но очень уж примитивная.
И как я ни крутил, ничего особо хорошего добиться не удалось. Вот это, пожалуй, лучшее, при том что понятно что обман:
(c = 0.2 - 1.428i, область 3x3).
Гипотеза: картинку "портит" полюс в нуле?
Zn+1 = Artcg(Zn)+ c?
Тоже не очень впечатляет (с = 1.15i):
Похоже, опять мешают полюса, которые у арктангенса сидят в ±i (нарисован квадрат 10 на 10). Но между этими полюсами, вероятно, должна же быть какая-то жизнь?
Последовательное увеличение подтверждает эту догадку:
(Центр квадрата смещён вниз на 0.35i)
Скажем, Zn+1 = Zni + c?
Как и ожидалось, влияние полюса в нуле оказывается довольно пагубным (с = -0.9i, квадрат 10 на 10):
Вертикальная линия -- следствие того, что операция возведения в степень на комплексной плоскости, вообще говоря, неоднозначна для нецелых степеней. Где-то должен проходить разрез по фазе. Библиотека комплексных вычислений, которую я написал 15 лет назад, почему-то провела его по π/2. Зачем -- сегодня уже фиг вспомнишь; в более поздние университские годы нас учили обычно проводить его вдоль нуля.
Злобный какой-то результат, вы не находите? :))